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I study how pulse to pulse phase coherence in a pulse train can survive super-broadening by 
extreme self phase modulation (SPM). Such pulse trains have been used in phase self-stabilizing 
schemes as an alternative to using a feedback process. However, such super-broadened pulses have 
undergone considerable distortion, and it is far from obvious that they necessarily retain any useful 
phase information. I propose measures of phase coherence applicable to such pulse trains, and use 
them to analyze numerical simulations comparable to self-stabilization experiments. 



I. INTRODUCTION 

In some recent experiments (e.g. (T], 0), carrier- 
envelope phase (CEP) stabilized idler pulses were gen- 
erated in an optical parametric amplifier (OPA) by com- 
bining a pump pulse with its super-broadened replica as 
the signal pulse. The phase relationships in the optical 
parametric interaction ensure that the phase of the idler 
pulse is independent of that of the pump. This princi- 
ple has been tested experimentally using a train of pump 
pulses with arbitrary CEP [J] , and a train of idler pulses 
with stabilized phase has been achieved. The distinguish- 
ing feature of this technique is that it acts on each pulse 
in isolation, and does not rely on feedback. The idler 
phase is stabilized to an unknown value determined by 
the details of the spectral broadening process. 

Fang et al. [l| broadened their signal pulses using ex- 
treme self-phase modulation (SPM) in sapphire, then 
selected suitable spectral components from the super- 
broadened spectrum to form the signal pulse for their 
OPA stage (see fig. [T]). The success of the experiment 
raises interesting questions about how the severely dis- 
torted temporal profile of the signal pulse can contain 
sufficient sensible phase information to achieve the de- 
sired result, especially given that small differences be- 
tween input pulses can be turned into large differences 
by the strong nonlinearity. 

In this paper I investigate how (inter-pulse) phase co- 
herence can survive in a train of pulses broadened by 
a third-order nonlinearity. In most pulse propagation 
models based on a complex envelope A{t), the nonlinear 
interaction term governing self-phase modulation (SPM) 
appears as x'-'^^l^l^^- This is linear in phase (of A), and 
is usually expected to be well-behaved, although intensity 
fluctuations in the pulse train lead to CEP fluctuations 
m the output However, the true form of the non- 

linear interaction is x^^^E^- In the envelope picture, this 
adds an extra (non resonant and CEP sensitive) term 
to the nonlinearity, i.e. x'^^^A'^. This initially leads to 
third harmonic generation (THG) which can disrupt the 
phase stabilization process; higher (odd) harmonics can 
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also be generated. Even more importantly, the THG is 
generated from the whole broadened spectrum, not just 
the initial central peak. Thus we cannot always assume 
that all of the THG contribution will remain safely out of 
range of our SPM broadened part: we might well expect 
that fluctuations will be strongly mixed into the nonlin- 
ear process, destroying any pulse-to-pulse repeatability. 
Since the true nonlinearity includes this wide range of 
multi-wave interactions, we need to be aware that all 
have the potential to obscure the phase coherence of the 
train of broadened pulses. 

I define criteria characterizing the effect of fluctuations 
in intensity, pulse width, and CEP on the phase of the 
broadened pulse. Applying these to simulations compa- 
rable to the experiment of Fang et al., I show the im- 
portance of minimising fluctuations in the input pulse 
train: even though SPM is insensitive to CEP variation, 
intensity and width fluctuations can play a significant 
role. The wide bandwidth of these pulses, along with the 
need to retain the full x^'^'' interaction, places significant 
demands on the required numerical resolution. There- 
fore simulations of pulse propagation were done using the 
PSSD method^], which evolve the carrier oscillations di- 
rectly. 

Because the strong nonlinearity means that the results 
will vary significantly from one set of parameters to an- 
other, the data presented in this paper serves largely 
to indicate the various issues involved in extreme SPM 
broadening, rather than being a set of specific predic- 
tions. Nevertheless, the simulation results indicate that 
the phase coherence in the broadened pulse train is main- 
tained only up to a point, and that sensitivity to pulse 
intensity or width fluctuations is much stronger than for 
CEP variation. Further, SPM can only broaden pulses 
by so much before multi-wave cross talk and/or THG 
destroy the pulse-to-pulse phase coherence. 

In section [H] the basic phase properties relevant for this 
analysis are described, and in section IIIII the measures 
of phase coherence are introduced. Section HVl describes 
the simulation methods used, and then in section IVl the 
effects of fluctuations in the input pulse train (intensity, 
pulse width, and CEP) are described. Finally, section IVTl 
comprises a brief discussion, and section lVlII contains my 
conclusions. 
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FIG. 1: Diagram of tlie Fang et al. experiment [ij. Pulses 
of unknown CEP from an input train have a small fraction 
split off at a beam splitter and then broadened in sapphire. 
The rest of the input pulse (i.e. the pump pulse) is then com- 
bined with the broadened portion, in a non-coUinear OPA, 
generating a phase stabilized idler pulse. 



II. PHASE PROPERTIES 

Here I focus on the spectral phase 4>s{uj), since a broad- 
ened pulse will suffer chirp and other distortion, remov- 
ing any straightforward way to characterize the phase in 
the time domain. This spectral phase will depend on the 
phase of the input pulse (j)p{uj), on its other properties 
(intensity, width, wavelength, etc), and on the propaga- 
tion mcdimn (e.g. sapphire). Considering just intensity 
/ and input CEP (pp variation, we have that the (broad- 
ened) signal pulse has a phase 

0,(/,0;a;) = 0p(L^) + $(/, 0; a;) (1) 
= 0p(w) -)-$(/o,0o;t^) 

+ AK/o,'^o)(/-/o)+^/H 

+ ~ (t>a) + 6^{uj), (2) 

where /q is the average intensity of pulses in the train, 
and 4>Q the average CEP; further terms including the ef- 
fect of pulse width (r) variation can easily be added. 
The function encodes the full effect of the complicated 
nonlinear interactions on the phase. Here the linear re- 
sponses of the phase about /q and (/)o are given by A/ and 
A^, while the remaining nonlinear response is contained 
in 5i, 5^. 

The broadened signal pulse is then mixed with the orig- 
inal pulse in an OPA. The difference- frequency idler pulse 
that results has a spectral phase which is the difference 
of the pump phase (f>p{uj) and the signal phase (j)s{uj), 

0j(/, 0;w) = (j)p{uj) - (t)s{l,(l)]uj) (3) 
= -$(/o,0o;w) 

-A/(/o,(/.o)(/-/o)-'5/(/,(/.;w) 
- A^(/o, (/.o)(0 - 0o) - 5c^{1. 4>\ ^)(-4) 

We can see that each idler pulse output from the OPA 
will have a phase (pi stabilized to a value set by the non- 
linear propagation through the sapphire ~ the unknown 
CEP 4>p{uj) of the input pump pulse has has canceled 
out. However, this phase (pi will vary from pulse to 
pulse depending on how the fluctuations in the input 
pulse train alter each instance of the linear and nonlin- 
ear evolution. Although 1 consider only the (three-wave) 
difference-frequency generating term, note that the full 



T^(2) OPA interaction contains potentially phase sensitive 
terms [gj]. 



III. PHASE COHERENCE 

In order to determine whether there is sufhcient useful 
pulse-to-pulse phase coherence after the SPM broaden- 
ing, I analyze the last four terms in eqn.(l2]) to define the 
necessary measures of phase coherence. 

I first consider the linear response of the output pulse 
to a variation in the input pulse; e.g. how the spectral 
phase of the pulse varies due to small changes in the in- 
tensity. Knowing the linear response enables us to set 
tolerances on the input pulse train to guarantee a cer- 
tain level of stability between pulses in the output train. 
We also need to know how reliable this estimated linear 
response is: since we are considering a complex nonlinear 
interaction, we cannot simply the response to input vari- 
ations will be linear. What we hope is that the dominant 
response will be linear, and that the nonlinear corrections 
will remain small. 

The quantities defined below in eqns. ([5]) and (O are 
expressed for simplicity only w.r.t. intensity fluctuations; 
but by swapping the intensity arguments /i, I2 for (e.g.) 
CEP's (ppi, (j)p2, we can equally well define a A^, cr^, ex- 
pressing the response to input CEP variation; just as by 
swapping the intensity Ii , I2 for pulse widths ri , T2 we 
could define A^, cr^. 

Since our simulations provide data (such as spectral 
phases) not easily accessible in experiments, we can use 
this to our advantage in determining the likely coherence 
properties of the broadened pulse train, and even suggest 
changes to improve performance. 



A. Linear response 

The linear response of the phase to intensity variation 
can be estimated by taking an ensemble average over 
a large number of simulations. Each simulation starts 
with a pulse selected at random from the distribution 
of all possible pulses in the input pulse train. I define 
the linear response as an ensemble average over pairs of 
broadened signal pulses: 



AKc^) = 



1/1-/2 1 



(5) 



More simply, we might calculate A/ by simply taking an 
average of phase differences over a suitable range of pulse 
variation. If this A/ is small, phase is little affected by 
input variation - at least to a linear approximation. For 
the case of intensity variation, this linear response has 
been measured experimentally, as reported by Li et al. 
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B. Nonlinear response 

I estimate the contribution of the nonhnear part of the 
phase response using 

\ [h-h) I 

Although <j'j{uj) is calculated in the same manner as a 
variance, it does not mean that (j)s varies randomly as the 
intensity (or some other parameter) shifts. This cr/(w) is 
simply a measure of how nonlinear the response of (ps is 
to pulse variation. 

If this cr|(w) is small (ct/ <C A/), then the phase 
changes due to pulse fluctuations is predominately lin- 
ear, and A/ is a useful quantity. Note that the size of 
CTj(tj) is strongly dependent on the range of variation in 
/; if that range is altered then cr|(w) should be recalcu- 
lated. 



C. Other measures of coherence 

The measures I introduce above could not be eas- 
ily measured in an experiment: although the spectral 
phase of pulses can be measured (e.g. fs*]), doing this for 
each individual pulse in a train of (dissimilar) broadened 
pulses would be a rather challenging task. Here they 
are intended primarily as a theoretical construct useful 
for analyzing simulation results. In contrast, Dudley and 
Coen [21 have proposed a measure of "shot to shot" coher- 
ence g^i2 J which can be calculated by taking an ensemble 
average of pairs of results taken from a set of simulations. 
The measure is 



{El{\,h)E2{\M)) 



\E,{\h)\ ) {\E2{\t2) 



7) 



Although this is easy to calculate from a simulation, 
it includes contributions from intensity variation in the 
broadened pulses, a complication which we need to avoid 
in this work. 



IV. SIMULATIONS 

The simulations were done using the PSSD method 5, 
IJ^, which offers significant advantages over the tradi- 
tional FDTD and Pseudospectral Time-Domain (PSTD) 
techniques (see e.g. [ll|) for modeling the propagation 
and interaction of few-cycle pulses. Run times are gener- 
ally faster, and the PSSD method also offers far greater 
flexibility in the handling of dispersion. Whereas FDTD 
& PSTD propagate fields E{z), H{z) forward in time, 
PSSD propagates fields E{t), H{t) forward in space. Un- 
der PSSD, the entire time-history (and therefore fre- 
quency content) of the pulse is known at any point in 



space, so arbitrary dispersion incurs no extra compu- 
tational penalty. In contrast, the FDTD or PSTD ap- 
proaches must use convolutions or time-response models 
for dispersion. 

However, since z-propagated simulations do not han- 



.(6 



die either reflections or backward waves easily, we need 



to ensure we remain in the uni-directional propagation 
limit, where for an n-th order perturbative nonlinearity, 
„^(n)£;"-i/n,2 ^ I ^ggg g g_ j^-)^ ]v^Q|.g ^Yso that it is 

also possible to do these simulations using explicitly di- 
rectional flelds [l^, [13, [31 or even wideband envelopes 

MM- 

I apply the PSSD algorithm to the source-free 
Maxwell's equations in non-magnetic media, with an in- 
stantaneous x^^'' nonlinearity, so that the equations for 
E and H in the ID (plane wave) limit are 



dHy{t; z) 

dz 
dExjt; z) 

dz 



eosrit) * E,{t; z) + eo/'^^E^it; z)^ 8,) 
dt I ^ 

d 

-— [nQHy(t-z)\ , 



(9) 



where the * denotes the convolution necessary to allow 
for the dispersion of the medium. 

Typical array sizes used were N = 2^^, with time res- 
olution of O.lfs (hence a time window of T = 3.2768ps). 
Spatial propagation steps were chosen to be dz = 
OAcT/N w 12nm in order to ensure numerical stability. 
The pulse proflle used as an initial condition was 



E{t) — Eq am{Lj pt + (j)) sech{t/T). 



(10) 



Pulses were propagated through 2mm of sapphire mod- 
eled u sing the dispersion parameters from DeFranzo & 
Pazol [171, and nonlinearity data from Major et al. [T§ |. 
The pulse wavelength was Ap = 786nm, which corre- 
sponds to u!p ~ 2.4 X lO^^rad/s. The reference peak 
intensity Iref for the part of the pulse being SPM broad- 
ened was chosen to be compatible with the pulse energy 
of 4% of 1.5mJ (length IBOfs) reported by Fang et al. 0], 
being hef = 0.33 x lO^W/cm^. In [l|, their OPA stage 
selected a a signal wavelength of As ~ 1400nm by angle 
tuning, generating an idler Ai 1600nm. Here, however, 
we simulate only the important SPM broadening stage. 

It is worth noting that in these simulations, a 
transform-limited 130fs pulse is seen to be too narrow- 
band to undergo sufflent broadening before multi-wave 
cross-talk destroyed the phase coherence. Consequently, 
we shortened our default pulse length to 30fs, as the sim- 
pler alternative to adding a strong chirp. 

The simulations did not incorporate transverse effects 
in order to keep computation times down. Each ID simu- 
lation through 2mm of sapphire takes about one hour on 
a 2.4GHz PC, so individual simulations involving trans- 
verse effects would in principle be tractable over days or 
weeks. However, the results here depend on multiple sim- 
ulation sets: a small number of 3D simulations could not 
have addressed the concerns of this paper. 
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V. FLUCTUATIONS IN THE PULSE TRAIN 

Some previous work [l^, showed how pulses broad- 
ened in experiments Uke those of Fang et al.pj] might pre- 
serve useful phase information. However, these simula- 
tions used artificial nonlinear strengths and pulse widths 
to ensure the presence of a strong SPM lobe directly at 
the desired signal frequency in order to guarantee good 
performance. They also included the OPA stage, demon- 
strating that it could preserve the phase information sup- 
plied by the SPM lobe. These indicated a phase stability 
at the output of the OPA stage similar to that seen by 
Fang et al., even though the parameters used were not 
strictly comparable. 

I now present simulation results that match the pulse 
intensity and material dispersion more closely. Fang et 
al. were unlikely to have been operating in the regime 
where a strong SPM lobe sat at their chosen signal fre- 
quency, since that was far into the wings of the broad- 
ened spectra. Unfortunately, the match is not perfect 
since the experimental parameters quoted by Fang et al. 
are insufficient to fully characterize a simulation. Fang 
et al. used a near-degenerate OPA setup, so were look- 
ing to broaden their input pulses from 2.4rad/fs all the 
way down to 1.2rad/fs; thus we are most interested in 
the results in the nearby frequency range a; ~ 1.3 rad/fs 

It is therefore difficult to be sure whether a phase- 
stable train of broadened pulses will arise, since we 
need the wing of the broadened spectra to be stronger 
than the underlying background of complicated phase- 
sensitive nonlinear processes. When doing the large set 
of simulations representing the variability of the input 
pulse train, three different pulse properties were altered 
in turn: intensity, pulse width, and CEP. We chose varia- 
tions of the order of a few percent for intensity and pulse 
width: specifically a range of 5% in steps of 0.5%. For 
CEP, we obtained results for a 7r/2 range in CEP in steps 
of 7r/16, which provides results applicable for the full 2tt. 

The spectral intensities at our chosen input intensities 
and propagation distances are shown on fig. [2l We see 
that the lowest input intensity pulse has only minimal 
power present in the range of interest, even for the full 
2mm propagation, whereas the highest intensity pulse 
broadens rapidly. Further, we see (in the high inten- 
sity case) the usual series of SPM lobes pushing down 
to lower frequencies as the pulse propagates. We need 
to bear these spectra in mind when analyzing the follow- 
ing results, since a well behaved phase response will be 
of little practical use if those frequency components are 
undetectable. This caveat is particularly relevant in the 
lowest pulse intensity cases. 

Note that in all the following figures, the central fre- 
quency of the input pulse is off the right hand side of the 
frame, since LOp ~ 2.4rad/fs. 




frequency a; (rad/fs) 




J I I , I I I I I L 



1.5 2 
frequency a; (rod/fs) 

FIG. 2: Spectral intensities for (a) Jo — Iref ~ 0.33 x 
lO^^'W/cm^, (b) lo = 3/re/; at (solid line, red, medium) 
0.85mm, (dashed line, green, light) 1.43mm, (dot-dashed line, 
blue, dark) 2.00mm. Peak intensity values ai z — 0mm are 
about 10^. 



A. Intensity fluctuations 

Pure SPM causes a phase shift proportional to inten- 
sity, but here I am instead interested in the unavoidable 
and potentially significant complications arising because 
the nonlinearity also includes THG effects. 

Fig. [3] plots the logarithm of the phase variance cr| 
cause by intensity fluctuations, for three propagation dis- 
tances at two different intensities. Linear phase response 
only holds to a level of O.lrad/% input intensity fluctu- 
ations where the log- variance is less than —2. In these 
results, it is important to note that if the intensity fluctu- 
ates, the central frequencies of the SPM lobes also shift, 
and this can cause large changes in the phase variance 

2 

aj. 

First consider the results for the reference intensity 
Iq = Iref, fig- ISIa). By comparing the logj^p '^i curves at 
different distances, we see that the trend at higher fre- 
quencies is for the phase variances to decrease, a result of 
the SPM broadened spectral peak pushing outwards. In 
contrast, for lower frequencies the opposite trend occurs, 
because the gradually increasing background of phase- 
sensitive processes still tends to dominate the spectral 
wings of the pulse. Indeed, for the 2mm propagation dis- 
tance, the variances have increased to the point where 
any hope of a linear phase response has been lost. Note 
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frequency cj (rad/fs) 




frequency cj (rad/fs) 



FIG. 3: The effect of intensity fluctuations on the phase 
coherence measure for (a) Iq = Iref = 0.33 X fO"W/cm^ (c) 
lo = 3/re/; (solid line, red, medium) 0.85mm, (dashed line, 
green, light) 1.43mm, (dot-dashed line, blue, dark) 2.00mm. 
A value of log(cr|) — —2 corresponds to a nonlinear phase 
adjustment 5i ~ ±0.1 rad / %. Width r = 30fs; Ap = 786nm; 
^ tjp ~ 2.4rad/fs. 



— 1 I I I I I I I r 

0.2 - 



<" - 



— 2 - '' -"^^iT^nn 



1.2 1.4 1.6 1.8 

frequency co (rad/fs) 

FIG. 4: The linear response of the phase to intensity fluctua- 
tions at 1.43mm and (dot-dashed blue line, dark) Iq = 3/re/. 
Width r = 30fs; \p = 786nm; ^ t.;^ ~ 2.4rad/fs. 

correspond to bandwidth fluctuations, they also cause 
different amounts of SPM-induced broadening, and shift- 
ing SPM lobe positions; hence we see similar trends as 
for intensity fluctuations. 

Again we see the spike at Iq — Iref, 1.3rad/fs and 
1.43mm, caused by the dip in the intensity spectra. Also, 
the variances tend to increase with propagation distance. 
Interestingly, the Jq = 3/re/ variances are smaller than 
those for the lower intensity Iq = Iref case, giving an ex- 
ample of how the balance between phase-sensitive nonlin- 
ear effects and coherent (SPM) spectral broadening can 
change. More generally, however, variations of ~ 1%) 
in width or intensity have a comparable effect, although 
there are many differences of detail. 



the spike in aj for 1.3rad/fs at 1.43mm, this corresponds 
to the dip in intensity seen on fig. [Ifa). 

Similar trends are shown for the higher intensity /q = 
3/re/, fig- [3Ic)j although the centre of the spectrum 
rapidly loses phase coherence, largely due to the strong 
SPM-induced spectral modulation. Note also the signif- 
icant modulation of the variances caused by the SPM 
lobes in the intensity spectrum, which are clearly visible 
early in the propagation at 0.85mm. 

For both intensities, there is a window of low phase 
variance at 1.43mm, although only for the /q = 3/re/ 
case is there also appreciable spectral intensity in this 
region. Therefore on fig. 2] we can plot A/ this (3/re/) 
case. For an ujg — 1.3rad/fs, we see that A/ ~ —0.1 
rad/%, hence phase stability is maintained to less than 
O.lrad if the intensity variation between the train of input 
pulses is less than 1%. However, in this particular case, 
the linear phase shift (for a 1% variation) is of the same 
order as the nonlinear contribution. 



B. Width fluctuations 

Fig. [5] shows the effect of pulse width fluctuation on 
the output phase response. Since the width fluctuations 



C. CEP fluctuations 

The final type of pulse variation we consider is CEP 
fluctuations. This case differs markedly from those for 
variations in intensity or width, since those were dom- 
inated by the concomitant alterations to the SPM. In 
contrast, since SPM is CEP insensitive, CEP fluctuations 
can only act (at least initially) through the THG-like con- 
tribution to the nonlinearity. 

The results in fig. [Slshow that even for the high inten- 
sity (/q = 3/re/) case, the response of the output spec- 
tral phase to the full range of CEP variation in the in- 
put pulse is remarkably linear. This emphasizes that the 
CEP stabilization scheme can work as expected, produc- 
ing a (nearly) input-CEP independent result - as long as 
the intensity and other parameters are sufficiently well 
stabilized. 

It is worth noting that although these results show 
that CEP effects are negligible in the regimes consid- 
ered in this paper, we still see that the variances tend 
to increase for either longer propagation distances or for 
higher pulse intensities. Indeed, for the high intensity 
case at z = 2.00mm, CEP fluctuations increase the log- 
variance to about —3, and have some non- negligible ef- 
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frequency cj (rad/fs) 




frequency co (rad/fs) 



FIG. 5: The efltect of pulse width fluctuations on the phase 
coherence measure, for (a) Jo = Iref = 0.33X lO'^'W/cm^ (c) 
Jo = 3/re/; (solid line, red, medium) 0.85mm, (dashed line, 
green, light) 1.43mm, (dot-dashed line, blue, dark) 2.00mm. 
A value of log(cr?)) — —2 corresponds to a nonlinear phase 
adjustment St ~ ±0.1 rad / %. Average width r = 30fs; 
Ap = 786nm; — > tjp ~ 2.4rad/fs. 

feet on the linear response to CEP fluctuations. Clearly, 
even a moderate further increase in intensity or propaga- 
tion distance would move the propagation into a regime 
where the train of broadened pulses is no longer coherent, 
so that the Fang et al. scheme could no longer generate 
phase stabilized pulse trains. 

VI. DISCUSSION 

This analysis has characterized the significance of CEP 
sensitive processes in the pulse broadening process. In 
particular, tr^ characterises the strength of these pro- 
cesses, and highlights important features of extreme- 
SPM propagation. However, since the usual expectation 
is that the CEP sensitive THG terms have minimal effect 
and the SPM term is dominant, it is worth considering 
why the results here do not always return an unambigu- 
ously linear phase response. 

Firstly, since the desired signal frequency will be in the 
wings of the broadened spectrum, small nonlinear effects 
can easily be significant. Secondly, the THG term ap- 
plies to the entire spectrum, including its low frequency 
wing. Both these contributions are enhanced by long 
propagation distances, which provide more than enough 




frequency cj (rad/fs) 




J I I I I I I I , L 



1 .5 2 

frequency cj (rad/fs) 

FIG. 6: The effect of input phase fluctuations on the phase 
coherence measure, for (a) 7o = Iref = 0.33xl0'''W/cm^ (c) 
Jo = 3/re/; (solid line, red, medium) 0.85mm, (dashed line, 
green, light) 1.43mm, (dot-dashed line, blue, dark) 2.00mm. 
A value of log(a"^)) = —2 corresponds to a nonlinear phase 
adjustment ~ ±0.1 rad / %. Width r = 30fs; Ap = 786nm; 
^ tjp ~ 2.4rad/fs. 

opportunity for CEP sensitive effects to accumulate, and 
then fold themselves back in to the propagating pulse. 

One might wish to compare simulations using a full 
T^(3)£;3 approach against an SPM-only model. Although 
non-trivial in standard PSSD, it is possible to use wide- 
band envelope techniques 15, 16] which allow the 
X^^^E^ term to be spht efficiently into SPM and THG 
parts. However, this introduces a non-physical disper- 
sion that averages over carrier cycles, which removes the 
apparent value of such comparisons. A more physical so- 
lution is to alter the dispersion above (e.g.) to ~ 2tt'o, to 
guarantee that there is no significant THG from the bulk 
of the pulse; however this leaves in place CEP sensitive 
THG from the low-frequency wing of the pulse spectra. 
Consquently it is not clear how to compare SPM and 
THG effects in a way that makes physical sense. 

VII. CONCLUSIONS 

In this paper I have proposed measures to assess the 
phase coherence of a pulse train subject to propagation 
through a nonlinear dispersive medium. I then used these 
to numerically investigate how pulse trains broadened by 
extreme SPM could retain pulse-to-pulse phase coher- 
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ence, as required by the Fang et al. [l| scheme. The 
A/ measure allows us to set tolerances in intensity fluc- 
tuations that guarantee a level of output stability, and 
the aj measure tells us how reliable those tolerances will 
be. The response to CEP {4> — > A^, cr^) and pulse width 
(t AtjCt^) fluctuations was also investigated. Note 
that such trains of broadened pulses retain remarkably 
good phase coherence in response to CEP fluctuations, 
but that there still remains an underlying sensitivity to 
CEP changes. Pulse with a wider bandwidth give a more 
robust spectral phase, since extra bandwidth allows for 
rapid broadening to the desired spectral range before 
multi-wave cross-talk can degrade (or destroys) the phase 
coherence. 

Simulations of the type used in this paper along with 
the (T^ measures can provide useful insight into experi- 
mental design. However, the simulation results presented 



here contain a great deal of fine detail, so such simula- 
tions need to be carefully customized to the desired ex- 
perimental parameters. This also suggests the potential 
for obtaining improved results in experiment by making 
rather small (but very specific) changes to the operating 
parameters. 

Lastly, the dominant process causing the loss of phase 
coherence in the broadened pulse train is intensity fluc- 
tuations in the input; varying the CEP of each pulse in 
the input train generally has a much smaller effect. 
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Appendix: CLEO'06 summary, "Phase Retention in 
SPM Super-broadened Pulses" 

P. Kinsler, S.B. Radnor, J.C.A. Tyrrell, G.H.C. New, 

Pulses that have been super-broadened by a third-order 
nonlinearity are frequently used when generating a phase 
stabilized output. But how can such temporally mangled 
pulses retain any useful phase information? 

We analyse recent experiments (such as [EIll) which 
phase-stabilized idler pulses are generated in an opti- 
cal parametric amplifier (OPA) by combining a pump 
pulse with its super-broadened replica. The success of 
this self-stabilization experiment raises interesting ques- 
tions about how the severely distorted temporal profile of 
the super-broadened signal can contain sufhcient sensible 
phase information to achieve the desired result. Our re- 
sults demonstrate the mechanism by which phase order 
survives amid the apparent temporal chaos of the pulse's 
electric field £^p„i,e- Our analysis is informed by nu- 
merical simulations using our recently-developed Pseudo- 
Spectral Space Domain (PSSD) technique @, and some 
results from one instructive example are shown in fig. [T] 
In traditional pulse propagation models based on the 
complex envelope A{t), the nonlinear interaction term 
governing SPM appears as x^^^'I^P^- This is linear in 
the phase, and might therefore be expected to be rela- 
tively well-behaved. However, one cannot model super- 
broadened pulses using envelopes, even with the various 
corrections to envelope propagation that are available. 
One must therefore resort to a model capable of evolv- 
ing the carrier oscillations directly, either the well known 
FDTD, or other methods like PSSD[| or G± variables 
[ist . Given that the form of the nonlinear interaction is 
now x^^'^E^, where E{t) is the complete field, it becomes 
unclear how the multiplicity of potential three-wave in- 
teractions could generate and retain the necessary phase 
information. The crucial questions we address are: 

• How does phase information survive in the presence 
of severe self-phase modulation and self-steepening? 

• Given that it does survive, how sensitive is the process 
to intensity and phase fluctuations in the pump pulse? 

• How does the number of cycles in the pulse affect the 
behaviour? 

Our simulations show firstly that a theoretical model 
of the SPM-OPA system described in [U, Q can repro- 
duce the character of the observed phase stability. Close 
analysis of the results shows that, despite the mangled 
appearance of the super-broadened pulse, some parts of 
it retain an ordered (phase) structure, and it is these sec- 
tions that are used in the experiment (e.g. the 1st lobe 
in fig. ED. 

The investigation has also revealed more general rules 
about the retention of phase information in strong SPM. 
The simulations have also enabled us to quantify the sen- 
sitivity of the phase-stabilisation process to fluctuations 
in the pump intensity, a feature that bears directly on the 
viability of the technique in the laboratory. The pump 
laser used in P, Q has an intensity stability to within 
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FIG. 7: Pulses at 786nm super-broadened in a sapphire-like 
medium to generate a lowest-frequency (1st) spectral lobe at 
~ 1200nm to be used in the OPA stage of a phase stabilization 
experiment il,, 2j]. • TOP; Starting at the bottom, this frame 
shows the pulse's electric field (iSpuiBc), a time-domain recon- 
struction of the electric field using just only the frequency 
content of this 1st lobe (i.e. Eli), and a similar reconstruc- 
tion {El2) from just the second lobe. The upper curves are 
the phase structure (j)Li and 4>l2 first and second spectral 
lobes respectively. • BOTTOM: Spectra of 30fs (top) and 
6fs (bottom) pulses at 786nm after similar amounts of spec- 
tral broadening, which require different propagation lengths. 
The nearly identical full and dotted lines compare the spec- 
tra differing by the influence of a 5% intensity variation. The 
other (offset) full lines shows the difference between those very 
similar pairs of spectra at the two intensities. 



1%, leading to absolute phase fluctuations of ^ O.Itt ra- 
dians. This compares remarkably well with our predicted 
0.237rad (see fig. [5]), especially given that our model pa- 
rameters do not closely match the experiment. Of course 
this match occurs because the phase stability (w.r.t. in- 
tensity) is still linear, and proportional to the broadening 
of the spectral component phase-matched to the signal 
frequency selected by the set-up of the OPA stage. Since 
these features are common to both our model and the ex- 
periment, the other differences of detail do little to upset 
the comparison. 

Lastly, we investigate the phase stability in the limit 
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of very few cycle pulses. There are two effects here that 
can shift the phase variation out of the linear regime seen 
for many cycle pulses. The first is the nonlinear response 
time, as modelled in [2l[, who used x^^'' with a Ifs re- 
sponse time; however, note that recent work [l8| gives a 
faster response (and thus will lead to smaller phase cor- 
rections). The second is the strong non-SPM effects from 
the x^^^ interaction. We show how these two nonlinear 
corrections to the phase properties affect the phase sta- 
bility in the few-cycle limit, and thus how they might 
affect experiments such as [l], Q if pushed to those ex- 
treme limits. 



FIG. 8: Phase stability as a function of percentage inten- 
sity variation for an initial 30fs pulse, as per fig. [7] The 
4>SPM data is extracted from the 1st spectral lobe after the 
SPM broadening process. The 4>idier data uses the original 
pump pulse and the SPM broadened copy in an OPA stage 
to generate a phase-stabilized idler. 
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